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Abstract. We discuss physical properties and baryonic content of broad Lya absorbers (BLAs) at low redshift. 
Fhese absorption systems, recently discovered in high-resolution, high-signal to noise quasar absorption line 
spectra, possibly trace the warm-hot intergalactic medium (WHIM) in the temperature range between 10 and 
10 6 K. The central idea is that in ionization equilibrium WHIM filaments should contain a very small fraction 
of neutral gas (/hi ~ 10 -5 — 10 -6 , typically), giving rise to weak intervening Hi Lya absorption. Due to the 
' high temperature of the WHIM, these Lya absorbers must be thermally broadened to Doppler parameters (6 

(f) , values) > 40 kms~ . It is expected, however, that also non-thermal line broadening processes, line blends, and 

■ noise features can mimic broad spectral features, complicating the quantitative estimate of the baryon content 

of the BLAs. To extend previous BLA measurements we have reanalyzed archival STIS data of the two quasars 
H 1821+643 and PG 0953+415 and have identified 13 BLA candidates along a total (unblocked) redshift path of 
Az — 0.440. Combining our measurements with previous results for the lines of sight toward PG 1259+593 and 
PG 1116+215, the resulting new BLA sample consists of 20 reliably detected systems as well as 29 additional 
tentative cases, implying a BLA number density of cJJVbla /<iz = 22 — 53. Eight BLAs show associated absorption 
|' from Hi Ly/3 and/or Ovi. The comparison between Lya, Ly/3, and Ovi line widths suggests that non-thermal 

1~ \ broadening and noise features substantially affect the observed BLA b value distribution. However, it remains 

unclear whether BLAs and O VI absorbers trace the same gas phase in the WHIM filaments. We estimate that 
the contribution of BLAs to the baryon density at z = is f2b(BLA)> 0.0027 /170 _1 for absorbers with log 
[iV(cm _2 )/6(kms -1 )] > 11.3. This number indicates that WHIM broad Lya absorbers contain a substantial 
fraction of the baryons in the local Universe. 
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ion to trace the WHIM at temperatures of T ~ 3 x 10 5 
K in the FUV regime. Oxygen is a relatively abundant 
element and the two available O vi transitions (located 
at 1031.9 and 1037.6 A) have large oscillator strengths. 
A number of detections of intervening WHIM O vi ab- 
sorbers at z < 0.5 have been reported in the litera- 
ture (Tripp, Savage, & Jenkins 2000; Oegerle et al. 2000; 
Tripp & Savage 2000; Chen & Prochaska 2000; Savage et 
al. 2002; Richter et al. 2004, Sembach et al. 2004; Savage et 
al. 2005; Tripp et al. 2005; Danforth & Shull 2005). These 
measurements imply a number density of O vi absorbers 
per unit redshift of dNovi/dz m 17 ± 3 for equivalent 
widths W x > 30 mA (Danforth & Shull 2005). Assuming 
that 20 percent or less of the oxygen is present in the form 
of O vi (/ovi < 0.2) and further assuming a mean oxygen 
abundance of 0.1 solar, the measured number density of 
Ovi absorbers corresponds to a cosmological mass den- 
sity of Q b (0 vi)> 0.0022 h 70 _1 . For the interpretation of 
this value it has to be noted that O vi absorption mainly 
traces gas with temperatures around 3 x 10 5 K, but not the 
million-degree gas phase which probably contains the ma- 
jority of the baryons in the WHIM. Very recently, Savage 
et al. (2005) have reported the detection of Ne vin in an 
absorption system at z « 0.2 in the direction of the quasar 
HE 0226-4110. Ncvin traces gas at T ~ 7x 10 5 K (in col- 
lisional ionization equilibrium) and thus is possibly suited 
to complement the O vi measurements of the WHIM in a 
higher temperature regime. However, as the cosmic abun- 
dance of Nevm is relatively low, Nevin is not expected 
to be a particularly sensitive tracer of the WHIM at the 
signal-to-noise levels achievable with current UV spectro- 
graphs. This is supported by the non-detections of inter- 
vening Nevm in other high S/N STIS data (Richter et 
al. 2004). X-ray absorption measurements are very impor- 
tant for studying the WHIM but are currently limited in 
scope because of the small number of available background 
sources and the relatively low spectral resolution of cur- 
rent X-ray observatories (FWHM~ 500 to 1000 kms" 1 ). 
Weak WHIM X-ray absorption by Ovn and Ovm has 
been observed with the Chandra X-ray observatory to- 
wards the quasar H 1821+643 (Mathur, Weinberg, & Chen 
2003). More recently, Nicastro et al. (2005) have reported 
the detection of Ovn absorption in a WHIM filament 
in the Chandra spectrum of the flaring blazar Mrk421. 
In the future, more sensitive X-ray observatories such as 
the X-ray Evolving Universe Spectrometer (XEUS) and 
Constellation- X will be of great importance to study the 
properties of the warm-hot intergalactic gas over the en- 
tire temperature range in both absorption and emission. 

Next to high-ion absorption from oxygen and other 
metals, recent observations with STIS (Richter et al. 2004; 
Sembach et al. 2004) suggest that WHIM filaments can be 
detected in Ly a absorption of neutral hydrogen. Although 
the vast majority of the hydrogen in the WHIM is ionized 
(by collisional processes and UV radiation), a tiny frac- 
tion (/hi < 10 -5 , typically) of neutral hydrogen should 
be present. Depending on the total gas column density 
of a WHIM absorber and its temperature, weak Hi Lya 



absorption at column densities 12.5 < log iV(Hi)< 14.0 
may arise from WHIM filaments and could be used to 
trace the ionized hydrogen component. The Lya absorp- 
tion from WHIM filaments is expected to be very broad 
due to thermal line broadening, resulting in large Doppler 
parameters of b > 40 kms -1 . Such lines arc generally 
difficult to detect, as they are broad and shallow. High 
resolution, high S/N FUV spectra of QSOs with smooth 
background continua are required to successfully search 
for broad Lya absorption in the low-redshift WHIM. The 
Space Telescope Imaging Spectrograph (STIS) is the only 
instrument that has provided such data, but due to the 
instrumental limitations of space-based observatories, the 
number of QSO spectra adequate for searching for WHIM 
broad Lya absorption (in the following abbreviated as 
'BLA') is very limited. So far, two sight lines observed with 
STIS towards the quasars PG 1259+593 (z om = 0.478) 
and PG 1116+215 (z om = 0.176) have been carefully in- 
spected for the presence of BLAs, and a number of candi- 
dates have been identified (Richter et al. 2004; Sembach et 
al. 2004). These measurements imply a BLA number den- 
sity per unit redshift of dNshA/dz « 20 — 90 for Doppler 
parameters b > 40 kms^ 1 and restframe equivalent widths 
W > 30 mA. The large range for cL/Vbla/ dz partly is 
due to the uncertainty about defining reliable selection 
criteria for separating spurious cases from good broad 
Ly a candidates (see discussions in Richter et al. 2004 and 
Sembach et al. 2004). Transforming the number density 
c^Abla / dz into a cosmological baryonic mass density, one 
obtains f2 b (BLA)> 0.004 h 70 _1 . This limit is only about 
10 percent of the total baryonic mass density expected 
from the current cosmological models (see above), but is 
above the limit derived for the intervening O vi absorbers 
(ri;,(Ovi)> 0.0022 h 70 _1 ; see above). Clearly, it is of great 
importance to extend the current sample of BLA candi- 
dates using recently obtained and archival STIS data to 
more reliably deduce the baryonic content of the WHIM 
at low z. 

In this paper we discuss the physical properties of 
the BLAs in the low-redshift Universe and review the 
observational constraints that limit the detectability of 
these systems. Moreover, we reanalyze the high S/N STIS 
spectra of the two QSOs H 1821+643 (z om = 0.297) and 
PG 0953+415 (z cm = 0.239) to search for WHIM broad 
Lya absorption along these two sight lines (Sect. 3). We 
combine our data with the results from previous mea- 
surements for PG 1259+593 and PG 1166+215 and dis- 
cuss physical properties and baryon content of the BLAs 
in Sect. 4. A summary of our study is given in Sect. 5. 

2. Broad Lya absorption 

2.1. Motivation 

Lya lines arising in the WHIM will be thermally broad- 
ened, with b values exceeding ~ 40 kms -1 . Since the gas 
is ionized, only a very small fraction of the hydrogen will 
be in neutral form and the H I column density is expected 
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to be small. Lya absorption from the WHIM therefore 
should occur in the form of very broad, shallow absorp- 
tion lines. Their shape would differ significantly from the 
shape of ordinary Ly a forest lines, which are usually rela- 
tively narrow (b ~ 25 kms -1 ; see Richter et al. 2004). For 
pure thermal broadening and assuming a collisional ion- 
ization equilibrium, one can estimate the total hydrogen 
content of each absorber directly from the measured line 
width (b value) and the Hi column density (see Sect. 2.2). 
However, in addition to thermal line broadening, non- 
thermal broadening mechanisms exist, such as the Hubble 
flow or gas flows and turbulence within the cosmic web fil- 
aments. Also unresolved line blends can produce absorp- 
tion features that are broad and shallow. For the analysis 
of interstellar and intergalactic absorption lines it is there- 
fore common to assume that the measured b value of an 
absorption line is composed of a thermal component, b t h, 
and a non-thermal component, & non -th, so that : 



HI Lya absorption; log N(HI)=13.0 



b = Vbth 2 + k 



non- th ' 



(i) 



Current FUV absorption line data typically have lim- 
ited spectral resolution and moderate signal-to-noise ra- 
tios (S/N). This means that some of the apparently broad, 
shallow absorber candidates in these spectra may be pro- 
duced by noise features, blends of multiple narrow ab- 
sorption lines, and continuum undulations rather than by 
WHIM gas. To reliably estimate the baryon content of the 
WHIM from broad Lya absorption it is therefore neces- 
sary to assess the effects of non-thermal broadening pro- 
cesses and data quality issues. In the following sections we 
discuss these various aspects in more detail to identify dif- 
ficulties for the interpretation of broad spectral features. 

2.2. Thermal line broadening 

If the width of an absorption line is governed solely by 
thermal broadening (i.e., 6 n on-th = 0), the measured b 
value of the line serves as a direct measure for the tem- 
perature of the gas: 



2kT 



0.1291 



kms 



(2) 



where k is the Boltzmann constant, m is the particle mass, 
and A is the atomic weight of the element. We demon- 
strate the effect of thermal line broadening of the Lya 
absorption in Fig. 1 : we model the expected shape (Voigt 
profile) of an H I Ly a absorption line with a constant H I 
column density of log iV(Hl)= 13.0 for Doppler param- 
eters b = 20,40,80, and 150 kms -1 . From equation (2) 
follows that T « 606 2 for pure thermal broadening. The 
b values of 20,40,80, and 150 kms -1 given in Fig. 1 then 
correspond to temperatures of log T « 4.4, 5.0, 5.6, and 
6.1. Clearly, it is much easier to detect Lya absorption 
from gas with T « 10 4 K than Lya absorption from gas 
with T w 10 6 K. The search for broad Lya absorption 
from the WHIM at temperatures T — 10 5 - 10 6 K there- 
fore requires spectral data with good S/N and a smooth 
background continuum. 
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b=20 km s" 1 
log N(HII)=16.4 



b=40 km s" 1 
log N(HII)=17.8 



b=80 km s" 1 
log N(HII)=19.5 



b=150 km s" 1 
log N(HII)=19.8 



-500 +500 -500 +500 

Velocity [km s _1 ] 

Fig. 1. Theoretical Voigt line profiles for Hi Lya absorp- 
tion with a constant Hi column density of log iV(Hl)= 
13.0 and b values of 20,40,80, and 160 kms -1 are shown. 
For pure thermal broadening and in collisional ionization 
equilibrium, these lines would trace gas with ionized gas 
column densities of log iV(Hll)= 16.4, 17.8, 19.5, and 19.8, 
respectively. 



If collisional ionization equilibrium (CIE) applies, one 
can directly calculate the hydrogen ionization fraction log 
/ H =log ((H°+H+)/H°) w log (H+/H ) as a function of 
the gas temperature (Sutherland & Dopita 1993; Richter 
et al. 2004) : 



log/ H (T) w -13.9 + 5.4 logT- 0.33 (logT) 2 



(3) 



Together with the measured Hi column density, iV(Hi), 
one can determine the total hydrogen column density via 
N(R) = /h^V(Hi). Under these assumptions, the chosen 
sequence of b values in Fig. 1 not only is a sequence of 
increasing temperature of the absorbing gas, but also a 
sequence of increasing ionized hydrogen column density 
(log JV(Hll) w 16.4, 17.8, 19.5, and 19.8). The correlation 
between higher gas content and weaker H I absorption has 
important observational consequences if the lines are ther- 
mally broadened: the broader and shallower a WHIM Hi 
Lya line is, the more baryonic material it traces. This 
means that it will be particularly difficult to detect those 
WHIM Lya absorbers that have the largest baryon con- 
tent, as many of these features may fall below the de- 
tection limit. Note that for low gas densities (nn < 10~ 5 
cm -3 ) photoionization by the UV background contributes 
to the ionization of the WHIM (see Sect. 2.6). Therefore, 
the total hydrogen column density derived from equation 
(3) under the assumption of CIE has to be regarded as a 
lower limit for N(H) in BLAs. 



2.3. Non-thermal line broadening 

The network of WHIM filaments represents cosmic struc- 
tures at scales of a several Mpc (e.g., Cen & Ostriker 1999; 
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Fig. 2. Artificial spectra of Ly a absorption in the wavelength range between 1220 and 1240 A are shown. The upper 
panel shows the initial synthetic spectrum (without noise) with five Lya absorbers A, B, C, D, and E having log 
N{Rl)= 13.23, 12.28, 13.07, 13.10, 13.12 and b values of 141,32,24,25,51 kms -1 , respectively. The lower three panels 
show the same spectrum with Poisson noise added at S/N levels of 50, 20 and 10 per 7 kms -1 wide pixel. For a 
S/N< 20 the broad Lya absorber A almost vanishes in the noise, and the two narrow absorbers C and D start to 
smear together to one broad spectral feature. At a S/N of 10, the broad Ly a absorber E is affected by a noise peak and 
thus looks asymmetric. These spectra demonstrate the difficulty to reliably identify intrinsically broad Lya absorbers 
in data with limited S/N. 



Dave et al. 2001), while the thickness of an individual fila- 
ment is expected to be on the order of several hundred kpc 
(Valageas, Schaeffer, & Silk 2003). Line broadening due 
to peculiar velocities and turbulence in the gas therefore 
must be considered when it comes to the interpretation 
of lines widths of WHIM absorbers. Clearly, at such large 
scales the gas is expected to have density and velocity 
substructure that will leave its imprint in the observed H I 
optical depth distribution (Kawahara et al. 2005). These 
effects would result in non-Gaussian shapes of broad H I 
lines, which could be identified in absorption line data 
at sufficient spectral resolution and S/N. Broad Hi lines 
that are asymmetric therefore should be interpreted with 
particular care, as their width may not represent a good 
measure of the temperature of the gas. 

As spectral resolution and signal-to- noise (S/N) in 
FUV data are limited, apparently broad Lya absorption 
may be caused by individual unresolved velocity com- 
ponents. In particular the measurements of line widths 
from intermediate and low spectral resolution spectra 
(FWHM> 20 kms -1 ) often turn out to be imprecise due 
to the presence of subcomponents that are not, or only 
barely resolved in the data. The resolution of the STIS 
E140M spectra used in this and the previous two studies 



is high, and with a FWHM of ~ 7 kms 1 , much smaller 
than the typical widths of the lines under investigation 
(b > 40 kms" 1 or FWHM> 66.5 kms" 1 ). At this resolu- 
tion, most of the existing subcomponent structure in inter- 
vening absorbers should be resolved. Multiple absorption 
components that are separated by <~ 20 — 40 kms -1 will 
in most cases produce an absorption structure that looks 
asymmetric and non-Gaussian. For uniformity, in this pa- 
per we concentrate on STIS E140M spectra, but we note 
that BLA candidates have also been identified in lower- 
resolution observations with the STIS first-order gratings 
(e.g., Bowen, Pettini, & Blades 2002). 

If the WHIM gas participates in the overall Hubble- 
flow, Hubble broadening will contribute with ~ 7 kms -1 
per hundred kpc thickness of the absorber to the total line 
width. However, if the thickness of the absorbing WHIM 
structures is on the order of a few hundred kpc (Valageas, 
Schaeffer, & Silk 2002), one expects that &Hubbie < 30 
kms -1 , typically. Moreover, as the WHIM is produced by 
collapsed gaseous structures that undergo gravitational 
shock-heating, one may argue that the resulting WHIM 
filaments primarily have decoupled from the Hubble flow 
during the collapse phase and that the gas dynamics is 
driven primarily by internal processes (this should be valid 
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at least for the main direction of the collapse). Hubble 
broadening therefore is not expected to dominate the 
widths of WHIM spectral lines. For those BLAs in which 
associated Ovi absorption is observed, one can evaluate 
the possible presence of non-thermal line broadening by 
comparing the b value of Hi and Ovi. For pure thermal 
broadening the Ovi b value is expected to be four times 
smaller than that of Hi owing to the 16 times larger mass 
of oxygen (see equation 2). Non-thermal contributions to b 
may be significant for BLA/Ovi absorber pairs that have 
&hi < 4 fro vi- However, only a small fraction of the BLAs 
show associated O vi absorption, so that this consistency 
check on the gas temperature can be applied only to a 
few cases. In addition, significant velocity offsets between 
H i and O vi as well possible multi-component structure 
in the Ovi are observed in some cases. It thus remains 
unclear whether the Hi and the Ovi absorption always 
trace the same gas phase, and whether the comparison 
between O VI and H I line widths provides reliable infor- 
mation about non-thermal gas motions. 

2.4. Noise 

FUV spectral data of extragalactic background sources 
obtained with STIS typically have S/N< 30 per ~ 7 
kms -1 wide resolution element. The presence of noise 
clearly hampers the reliable detection of WHIM BLA can- 
didates. The S/N in the data will determine the detection 
limit of broad features and so it is clear that only the high- 
est H I column density WHIM filaments can be detected 
with current FUV data. In contrast to the O vi absorbers, 
the detection limit for broad Ly a absorption cannot eas- 
ily be quantified in terms of a minimum equivalent width, 
as both the depth and the width of a line limit its de- 
tectability. For instance, a relatively narrow Hi Lya line 
with b = 40 kms" 1 and log N = 13.0 (see Fig. 1, upper 
right panel) will have an equivalent width of ~ 50 mA and 
can easily be detected with STIS data at a S/N> 15. At 
this S/N, however, a very broad line with b = 150 kms -1 
and log N = 13.0 (Fig. 1, lower right panel) would not be 
detectable, although the equivalent width would be the 
same. From the analysis of STIS data we find that BLA 
candidates can be identified if the following condition is 
fulfilled: 

iV(HI) . 3 x 10 12 _ 2/1 /A . 

im-isw cm (kms } ■ (4) 

In this equation, S/N is the local signal-to- noise ratio per 
resolution element. The S/N ratio may vary significantly 
within a spectrum due to large-scale flux variations in the 
QSO spectrum and the wavelength-dependent sensitivity 
of the detector. The average S/N in our STIS data is <~ 15, 
so that the typical detection limit for detecting broad Ly a 
absorbers in these data is log (N/b) m i n > 11.3. 

Noise features will also influence the shape of spectral 
features and thus will complicate the search for WHIM 
BLA candidates. Therefore, we have studied the influence 



of noise features (Poisson noise) in a set of artificial spec- 
tra with various S/N levels. The FWHM in these spec- 
tra is similar to those in the STIS data and the pixel 
size corresponds to one resolution element (7 kms -1 ). 
An example spectrum is shown in Fig. 2. The outcome 
of our noise analysis of these mock spectra is the follow- 
ing: (1) noise can create structures in the continuum that 
could be mis-identified as intervening broad Ly a absorp- 
tion, (2) noise can create broad absorption features out 
of multiple, individual velocity components that are sep- 
arated only by ~ 20 — 50 kms -1 , (3) noise can substan- 
tially modify the shape of an existing broad absorber, so 
that Gaussian-shaped absorption lines become asymmet- 
ric, and vice versa. The binning of several pixels can help 
to check on the validity of an absorption feature, but we 
find that at S/N< 20 more than 50 percent of all apparent 
broad absorption features are caused by noise rather than 
by intrinisically broad lines. Next to Poisson noise, other 
systematic noise sources can be important for the inter- 
pretation of spectral features in FUV absorption line data. 
For the STIS instrument, possible order-edge artifacts and 
features caused by the obscuration by the field electrode 
of the MAMA detector (the "repeller wire" ; see Woodgate 
et al. 1998) have to be considered. The limited S/N in the 
currently existing FUV data introduces the largest uncer- 
tainty for the proper identification and analysis of BLA 
candidate lines. Note, however, that the broad lines chosen 
for this study currently represent the best ones available 
for the search of BLAs. 

2.5. Continuum undulations 

Due to the large width and the shallowness of the BLAs, 
one may argue that part of these features are not pro- 
duced by intervening absorption but rather arc caused 
caused by small-scale undulations in the overlaying quasar 
continuum. One simple way to test for the presence of 
continuum undulations in the STIS data is to analyze 
the smoothness of the quasar continuum shortwards and 
longwards of the QSO Lya emission. If some of the fea- 
tures identified as broad Lya absorbers are instead pro- 
duced by continuum undulations, one should expect that 
such undulations are distributed all over the spectrum, 
whereas true intergalactic BLAs can occur only in the 
wavelength range between AeLA.min = 1215.67 A and 
ABLA,max = (1 + ^qso) x 1215.67 A. However, in none 
of the analyzed spectra do broad features (if significant 
with respect to the ambient noise) occur at wavelengths 
ABLA.max > (1 + ^qso) x 1215.67 A. This implies that con- 
tinuum undulations do not interfere with our BLA search. 

Continuum undulations also may be caused by weak 
emission lines from the QSO. Possible candidate lines 
for the range A B LA,max < (1 + ^qso) x 1215.67 A are, 
for instance, Hi Ly/3A1025.7, OviA1031.9, CmA977.2. 
However, none of the BLA candidates considered in this 
paper lie close to the wavelengths at which QSO emission 
in these (and other) lines would be expected. 
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2.6. Deviations from CIE 

The determination of /h and iV(Hll) using equation (3) 
is valid only under the assumption that the gas is in 
collisional ionization equilibrium, which may not be the 
case in general. In fact, it is expected that for low- 
density gas (riH < 1CT 5 cm~ 3 ) photoionization from 
the UV background becomes important, so that the ion- 
ization fraction (as defined in Sect. 2.2) in low-density 
BLAs should be higher compared to what is estimated 
for pure collisional ionization. To investigate the influ- 
ence of photoionization we have analyzed the output of 
a WHIM simulation presented by Fang & Bryan (2001). 
This simulation includes a combined (time equilibrium) 
collisional ionization and photoionization model (hereafter 
referred to as CIE+PH), based on ionization fractions 
from Mazzotta et al. (1998) and the CLOUDY photoion- 
ization code (Ferland et al. 1998). The WHIM simulation 
predicts that for T = 10 5 - 10 7 K and n H = 10~ 6 - 10~ 4 
cm~ 3 the CIE+PH ionization fraction is (on average) ~ 50 
percent higher than for CIE alone. If we restrict the den- 
sity range to rin > 10~ 5 cm~ 3 , the CIE+PH ionization 
fraction typically lies ~ 15 percent above the CIE esti- 
mate. This shows that photoionization is important for the 
BLAs. However, since we do not have information about 
the actual hydrogen volume densities in these absorbers, 
we only can notice that with the CIE assumption we prob- 
ably underestimate the ionized hydrogen column density 
(and thus the baryon content) of BLAs by ~ 15 — 50 
percent. To more precisely pinpoint the role of photoion- 
ization in BLAs a more detailed analysis of broad Lya 
absorption in cosmological simulations is required. This 
will be presented in a future paper (Richter et al. 2005, in 
preparation) . 

2.7. The baryonic content of broad Lya absorbers 

The cosmological mass density of the WHIM broad Lya 
absorbers in terms of the current critical density p c can 
be estimated by 

fib (BLA) = fan N ( m h /E A ^< ( 5 ) 

Pc c ' 

with p, = 1.3, m H = 1.673 x 10~ 27 kg, H = 70 
kms" 1 Mpc -1 , and p c = 3H 2 /8ttG. The index i denotes 
an individual broad Lya system along a line of sight j. 
Each measured absorption system i is characterized by its 
neutral hydrogen column density, A(H l)y and ionization 
fraction, fn,ij (see equation 3). Each line of sight j has a 
characteristic redshift range Az = z max — z m i n in which 
broad Lya absorption may be detected. The comoving 
path length AXj available for the detection broad Lya 
systems then is given by: 

AX = 0.5 {[(l+;w-AV2) 2 -lH(l+z min +A 2B /2) 2 -l]}, (6) 

where we refer to a cosmology with q = 0. Due to the 
presence of Ly a forest lines and interstellar absorption 
lines a blocking correction Azb for each sight lines has to 
be applied. 



2.8. Summary: interpretation of Hi line widths 

As there currently are only very few techniques that can 
be used to constrain the baryon budget of the WHIM by 
direct observations, the analysis of broad Lya absorbers 
in low-z QSO absorption line spectra is of fundamental 
importance for our understanding of the distribution of 
baryonic matter in the local Universe. Every WHIM fila- 
ment that contains sufficient amounts of neutral hydrogen 
will produce a broad Lya absorption line, but, vice versa, 
not every broad feature observed in FUV data is caused by 
the WHIM. Moreover, the quantitative estimate of A(H n) 
from A (Hi) and 6(Hi) is hampered by the various alter- 
native broadening mechanisms discussed above. Therefore 
it is clear that the interpretation of broad absorption fea- 
tures in FUV spectra and the determination of Q& from 
BLA candidates is afflicted with systematic uncertainties. 
This should be kept in mind for the second part of the 
paper, where we analyze BLA candidates along four lines 
of sight and estimate their baryon content. 

3. STIS observations of low-z broad Lya 
absorbers 

3.1. Overview 

STIS data with medium- to high spectral resolution and 
good S/N are available only for a very limited number 
of QSO sight lines, as it requires a substantial amount 
of observing time to observe objects with V magni- 
tudes > 15 at high spectral resolution. In this study, 
we analyze broad Lya absorbers along the two lines of 
sight towards H 1821+643 and PG 0953+415 and com- 
bine these data with previously published measurements 
for PG 1259+593 (Richter et al. 2004) and PG 1116+215 
(Sembach et al. 2004). Information about the line-of-sight 
properties is summarized in Table 1. All four spectra were 
obtained between 1998 and 2002 using the medium res- 
olution FUV echelle mode (E140M), which provides a 
spectral resolution of A/AA « 46,000, corresponding to 
FWHMw 7 kms^ 1 . The integration time for each spec- 
trum varies between ~ 25 and ~ 95 ks (see Table 1). More 
detailed information about the observations and primary 
data reduction of these spectra can be found in Tripp, 
Savage, & Jenkins (2000), Tripp & Savage (2000), Richter 
et al. (2004), and Sembach et al. (2004, 2005a). Note that 
the STIS data sets for H 1821+643 and PG 0953+415 have 
been extensively used to study the distribution and prop- 
erties of intervening O vi absorbers at low redshifts (Tripp, 
Savage, & Jenkins 2000; Tripp & Savage 2000; Tripp et 
al. 2001; Savage et al. 2002), as well as to investigate the 
statistical and physical properties of the low— z Lya for- 
est (Dave & Tripp 2001). In the latter study, the authors 
compared the properties of the Lya forest in these two 
spectra with those of mock QSO spectra from hydrody- 
namical simulations. For the analysis, Dave & Tripp used 
an automated continuum fitting routine together with an 
automated Voigt profile fitting algorithm (AutoVP; see 
Dave et al. 1997) to efficiently handle the large number of 
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Table 1. Log of STIS observations 



Sight line 


I 


b 


V 




Grating 


Exposure time 


Broad Lya analysis 








[mag] 






[ks] 




H 1821+643 


94.00 


+27.42 


14.24 


0.297 


E140M 


25.5 


this study 


PG 0953+415 


179.79 


+51.71 


15.32 


0.239 


E140M 


24.5 


this study 


PG 1116+215 


223.36 


+68.21 


14.80 


0.177 


E140M/E230M 


35.5 


Sembach et al. (2004) 


PG 1259+593 


120.56 


+58.05 


15.84 


0.472 


E140M 


95.5 


Richter et al. (2004) 



lines in their data. With this method, they find 89 Lya 
forest absorbers with b < 40 kms" 1 along a total redshift 
path of Az = 0.335 (Dave & Tripp 2001; their Fig. 6), but 
apparently only two systems with b > 40 kms -1 . However, 
our experience with automated line fitting routines is that 
broad and shallow spectral features often are not identi- 
fied by these routines, or that they are inappropriately fit- 
ted, e.g., by multiple, narrow lines. A careful re-analysis 
of these data focussed on the search of broad and shallow 
absorption features is justified. 

3.2. Data analysis 

For the data analysis of the STIS spectra of H 1821+643 
and PG 0953+415 we follow the procedures from previous 
spectral anlyses of STIS data for PG 1259+593 (Richter 
et al. 2004) and PG 1116+215 (Sembach et al. 2004). From 
a careful by-eye inspection of the STIS data (put onto a 
heliocentric velocity scale) we first have created a cata- 
log of broad spectral features that represent candidates 
for intervening broad Hi Lya absorption (A = 1215.67 
A; oscillator strength / = 0.460; Morton 2003). As we 
do not want to include absorbers that possibly are asso- 
ciated with the background quasar itself, we have consid- 
ered only those lines that are more than 5000 kms" 1 away 
from the quasar. We have omitted from our candidate list 
those absorption features that exhibit clear evidence for 
sub-component structure, line blending, and irregularities 
due to continuum undulations. Moreover, we do not con- 
sider broad absorption features that are part of complex 
multicomponcnt absorption systems, as the width of these 
absorbers is particularly difficult to determine. We there- 
fore ignore hot gas that is situated in multiphase absorbers 
that may preferentially be related to gas in the vicinity of 
galaxies and galaxy groups. This introduces a selection 
bias in our BLA sample. Note that recent observations 
of Nevni and Ovi in a multiphase absorber at z w 0.2 
toward HE 0226—4110 suggest that such systems may con- 
tain significant amounts of hot gas (Savage et al. 2005). It 
is important to emphasize that our selection procedure 
is based solely on the visual inspection of absorption fea- 
tures in the spectrum. This method thus involves a certain 
degree of subjectiveness for the identification of BLA can- 
didates, in particular if it comes to the evaluation of noise 
features, line asymmetries, and continuum undulations. 
As outlined earlier, automated fitting procedures currently 
do not represent a reliable alternative to assess the signif- 
icance of broad spectral features. Therefore, we here rely 



on selection criteria based on our experience with previ- 
ous absorption line data. Note that the BLA candidates 
presented in this paper have been identified and evaluated 
independently by the various authors of this paper. 

After having selected a sample of BLA candidate lines 
we have cross-checked whether these features could be 
higher Lyman lines or metal lines belonging to other inter- 
galactic absorption systems, and modified the list accord- 
ingly. We then checked for associated broad Ly f3 and O VI 
absorption for the BLA candidates. BLAs that show asso- 
ciated Ly/3 absorption with 6Ly/3 < 40 kms" 1 have been 
excluded from our sample (see Sect. 4.2.1). For the final 
BLA candidate sample, the lines have been continuum- 
normalized. The continuum fitting was done for each line 
individually, using low-order polynomials that fit the con- 
tinuum around a line in range from approximately —1000 
to +1000 kms" 1 . The base points for the polynomial fit 
were obtained from a binned (10 pixels, typically) version 
of the spectrum to smooth out the noise in the continuum 
flux and thus to increase the reliability of the continuum 
fit. The normalized broad Lya candidate lines then were 
fitted with Voigt profiles convolved with a Gaussian line 
spread function that corresponds to the spectral resolution 
of the STIS instrument (FWHMw 7 kms" 1 ). As the lines 
under analysis generally have FWHM> 50 kms -1 , possi- 
ble small deviations of the STIS line spread function from 
a single-component Gaussian have no measurable influ- 
ence on the outcome of the profile fit. For the fitting pro- 
cedure we have used the Voigt profile/ x 2 minimization 
routine FITLYMAN installed in the ESO-MIDAS soft- 
ware package (Fontana & Ballester 1996), which deliv- 
ers line centroids, redshifts, equivalent widths, logarithmic 
column densities, b values, and their corresponding errors. 
The total uncertainty for N and b have been calculated by 
quadratically adding the error estimates from the fitting 
process to those from the continuum placement. Due to 
the fact the S/N varies over the spectrum, the detection 
limit for BLAs, log (N/b) m i n , as defined in Sect. 2.4, also 
varies across the spectrum. Many of the BLA candidates in 
our sample have a low S/N in their profile. In some cases, 
there is weak (but not compelling) evidence for asymme- 
tries and unresolved component structure. Although these 
systems are all clearly detected with W r /a > 4.4 except 
for one line with W T /a — 2.6, it is probable that some 
of them are not related to thermally broadened Hi but 
are caused by blends of several subcomponents and noise 
features. As the S/N in our data is not good enough to reli- 
ably evaluate the true component structure in these broad 
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Fig. 3. The well-detected broad Lya absorber at z = 0.18047 in the STIS spectrum of H 1821+643 is shown. This 
system has an Hi column density of log N(Hl)= 13.21 and a b value of ~ 56 kms" 1 (see also Table 2). The shape of 
this BLA differs significantly from the shape of the narrow z — 0.17924 Lya forest absorption near 1433.5 A. 



lines, we regard these absorbers as tentative BLA detec- 
tions. Throughout the following we label these cases ac- 
cordingly. Measured column densities, b values, and other 
quantities for each absorber in our final BLA candidate 
sample are listed in Tables 2 and 3. In Table 4 we list es- 
timates for the temperature, ionization fraction, and total 
hydrogen column density for each system (see Sect. 4.3). 
Statistical information for each sight line is presented in 
Table 5. 



3.3. H 1821+643 

The wavelength range available for detecting intervening 
broad Lya absorption towards H 1821+643 (z cm = 0.297) 
is from - 1230.0 A to - 1556.5 A (> 5000 kms" 1 away 
from the quasar), corresponding to a redshift interval of 
Az = 0.269. The typical S/N in this range is ~ 21 per 
resolution clement, so that the detection limit for broad 
Lya absorption (see equation 4) is log (N/b) ~ 11.2. 
The blocking of certain spectral regions due to other 
interstellar and intergalactic lines reduces the total un- 
blended redshift interval to a value of Az = 0.238. We find 
six reliably detected single-component BLA candidates 
at z = 0.02642,0.12221,0.16352,0.18047,0.22616, and 
0.26659. Equivalent widths, column densities, and b val- 
ues for these absorbers are listed in Table 2. The systems 
at z = 0.16352, 0.18047, and 0.22616 have been previ- 
ously identified as H I Ly a absorption in lower-resolution 
GHRS G140L data (Tripp, Lu, & Savage 1998). Two of 
these systems (at z = 0.22616 and 0.26659) show associ- 
ated Ovi absorption, which has been previously analyzed 
by Tripp, Savage & Jenkins (2000). For the z = 0.26659 
system we find associated Ly/3 absorption, which is also 
broad (b « 42 kms -1 ). Considering all six BLA candi- 
dates we obtain a BLA number density per unit redshift 
of dN BL /dz w 25 for log (N/b) > 11.2. As an example 
for broad Lya absorption towards H 1821+643 we show 
in Fig. 3 the spectral region between 1430 and 1440 A. 
The BLA near 1435 A (z = 0.18047) is clearly visible, and 



there is no evidence that this feature is produced by a con- 
tinuum undulation or a blend of several narrow absorption 
lines. Continuum normalized velocity profiles of all BLA 
candidates towards H 1821+643 are shown in Fig. 4. From 
Voigt profile fitting of the Lya profiles we derive b val- 
ues between 44 and 56 kms -1 and logarithmic Hi column 
densities between 13.19 and 13.65 (see Table 2 for details). 
All the BLAs candidate lines toward H 1821+643 therefore 
are relatively narrow, which may imply that some of them 
are not related to the WHIM, but rather are Lya forest 
lines broadened by non-thermal processes. Our measured 
values of b = 45.6+3.9 kms -1 and log N = 13.65+0.06 for 
the BLA at z = 0.26659 agree well with the numbers cited 
by Tripp, Savage, & Jenkins (2000) (6 = 44.6±gj kms" 1 
and log TV = 13.62 + 0.03) for this system. Note that we do 
not consider here the broad Ly a absorption found in the 
z = 0.12120 system, as this line is part of a complex mul- 
tiphase absorption system with very broad O vi absorp- 
tion, possibly indicating unresolved component structure 
and non-thermal broadening processes (Tripp et al. 2001). 
A complete IGM analysis of the H 1821+643 sight line will 
be provided in a separate paper (Sembach et al. 2005a). 



3.4. PG 0953+415 

Towards PG 0953+415 we detect seven BLA candidates at 
z = 0.02108, 0.04382, 0.09315, 0.10294, 0.17984, 0.19142, 
and 0.19360 in the available wavelength range between 
- 1220.0 and - 1485.9 A (Az = 0.219). After applying 
a blocking correction, the total redshift path for finding 
BLAs towards PG 0953+415 reduces to Az = 0.202. The 
systems at z = 0.02108, 0.10294, and 0.19142 are tentative 
detections (see Table 2, last column for a more detailed 
description of these lines). Another broad spectral feature 
is observed near 1282 A, corresponding to z = 0.05461 for 
Hi Lya. This feature most likely is an order-edge arti- 
fact caused by the slightly imperfect ripple correction in 
the STIS data (see Sect. 2.4) and will not be considered 
any further. The S/N in the PG 0953+415 spectrum is 
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Fig. 4. Broad Lya absorbers towards H 1821+643, plot- 
ted on a restframe velocity scale. Voigt-profiles fits are 
indicated with the solid line. Other absorption features 
are related to the Lya forest and metal lines. 
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~ 12, implying a typical detection limit for broad Lya 
absorption of log (N/b) ~ 11.4 (equation 4). The systems 
at z = 0.09315, 0.17984, 0.19142, and 0.19360 have been 
previously identified as Hi Lya by Savage et al. (2002). 
None of the broad Ly a candidates towards PG 0953+415 
shows associated O vi absorption; the two O vi absorbers 
identified by Savage et al. (2002) are both multiphase ab- 
sorbers with Hi Lyman series absorption having b < 40 
kms -1 . The system at z — 0.19360 is also detected in 
broad Ly/3 absorption (with a line width consistent with 
that of the Lya absorption). From the four reliably de- 
tected BLA candidates and the blocking-corrected red- 



Fig. 5. Broad Lya absorbers towards PG 0953+415, plot- 
ted on a restframe velocity scale. Voigt-profiles fits are in- 
dicated with the solid line. Other absorption features are 
related to the Lya forest and metal lines. 

shift path of Az = 0.202 we obtain dN BL /dz w 20 for log 
(N/b) > 11.3 for the line of sight towards PG 0953+415. 
Including the three tentative cases we find dNBL/dz « 35. 
All seven systems are shown in Fig. 5. Voigt profile fitting 
of these lines results in b values between 40 and 122 km s _1 
and logarithmic Hi column densities between 12.88 and 
13.91, as listed in Table 2. 
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Table 2. BLA candidates towards H 1821+643, PG 0953+415, PG 1116+215, and PG 1259+593 



Aobs 


z 


W r a 


log N a 


6 a 


log (N/b) 


Other det. 


Status b 




[A] 




[mA] 




[kins" 1 ] 




species 






H 1821+643 


1247.789 


0.02642 


79 + 10 


13.24 + 0.07 


51.3 + 6.0 


11.5 








1364.232 


0.12221 


74 ±8 


13.23 + 0.05 


44.4 ± 3.8 


11.6 








1414.461 


0.16352 


75 ±9 


13.19 + 0.06 


55.7 + 4.9 


11.4 








1435.057 


0.18047 


68 + 9 


13.21 + 0.08 


52.8 + 5.4 


11.5 








1490.605 


0.22616 


139 ± 14 


13.47 + 0.07 


53.8 + 5.7 


11.7 


(Ovi) 






1539.750 


0.26659 


174 ± 13 


13.65 + 0.06 


45.6 ± 3.9 


12.0 


Ly/?, Ovi 






PG 0953+415 


1241.302 


0.02108 


69+11 


13.12 + 0.10 


65.1 + 9.1 


11.3 




tentative, S 




1268.938 


0.04382 


119 + 15 


13.34 + 0.07 


47.2 + 4.2 


11.7 








1328.905 


0.09315 


194 ± 16 


13.66 + 0.06 


43.9 ± 5.3 


12.0 








1340.807 


0.10294 


114 + 18 


13.32 + 0.16 


121.7 + 22.4 


11.2 




tentative, S 




1434.275 


0.17982 


74 + 10 


13.24 + 0.07 


46.0 + 5.5 


11.6 








1448.378 


0.19142 


126 + 10 


13.47 ± 0.08 


43.8 + 3.7 


11.8 




tentative, A 




1451.026 


0.19360 


232 ± 10 


13.91 ± 0.04 


40.3 + 2.8 


12.3 


Ly/? 






PG 1116+215 


1235.550 


0.01635 


111 ± 10 


13.39 + 0.06 


48.5 + 5.1 


11.7 








1265.820 


0.04125 


78 + 16 


13.24±g;ff 


105.0 + 18.0 


11.2 


Ovi 


tentative, M 




1289.490 


0.06072 


80 + 8 


13.28 + 0.06 


55.4 + 5.8 


11.5 








1291.580 


0.06244 


74 + 9 


13.18 + 0.07 


77.3 + 9.0 


11.3 


Ovi 






1320.060 


0.08587 


18 + 7 


12.90 + 0.10 


52.0 ± 14.0 


11.2 




tentative, S 




1328.470 


0.09279 


111 ± 14 


13.39 + 0.09 


133.0 ± 17.0 


11.3 




tentative, B 




1378.210 


0.13370 


86+11 


13.27 + 0.08 


83.6 + 10.4 


11.3 




tentative, A 




PG 1259+593 


1218.449 


0.00229 


190 ± 24 


13.57 + 0.10 


42.1+4.4 


11.9 


Ly/?, Ovi 






1221.024 


0.00440 


215 + 24 


13.65 + 0.14 


142.6 ± 29.4 


11.5 




tentative, S 




1221.812 


0.00505 


133 ± 20 


13.44 + 0.10 


57.2 + 8.9 


11.7 




tentative, S 




1232.046 


0.01347 


90+17 


13.24 + 0.14 


104.1 ± 17.5 


11.2 




tentative, S 




1267.334 


0.04250 


58 + 7 


13.06 + 0.06 


40.7 + 4.9 


11.5 




tentative, A 




1283.572 


0.05586 


149 ± 20 


13.46 + 0.14 


144.2 ± 27.9 


11.3 




tentative, A 




1298.456 


0.06810 


72 + 10 


13.14 + 0.09 


87.7 + 9.2 


11.2 




tentative, M 




1313.423 


0.08041 


45 + 7 


12.97 + 0.10 


42.0 + 4.5 


11.3 








1327.468 


0.09196 


72 ± 22 


13.13 + 0.21 


113.4 ± 29.2 


11.1 




tentative, A 




1340.654 


0.10281 


136+18 


13.41 + 0.17 


196.9 + 31.2 


11.1 




tentative, A 




1377.977 


0.13351 


53+11 


13.01 + 0.08 


48.1+4.3 


11.3 




tentative, A 




1386.276 


0.14034 


60+10 


13.06 + 0.07 


56.7 + 5.4 


11.3 




tentative, S 




1390.498 


0.14381 


80+15 


13.18 + 0.13 


115.8 + 9.6 


11.1 




tentative, S 




1396.225 


0.14852 


254 + 9 


13.91 + 0.06 


42.1 + 2.4 


12.3 


Ly/? 






1399.668 


0.15136 


106 + 9 


13.32 + 0.09 


65.3 + 5.5 


11.5 








1418.044 


0.16647 


113 + 13 


13.34 + 0.08 


93.3 + 8.9 


11.4 




tentative, A, 


R 


1424.130 


0.17148 


94 + 19 


13.25 + 0.16 


131.5 ± 17.4 


11.1 




tentative, M, 


R 


1433.167 


0.17891 


101 ± 10 


13.29 + 0.10 


98.5 + 9.1 


11.3 




tentative, M 




1441.174 


0.18550 


77+11 


13.17 + 0.12 


86.4+10.0 


11.2 




tentative, A, 


R 


1493.587 


0.22861 


133 + 10 


13.47 + 0.05 


40.3 + 2.9 


11.9 








1508.963 


0.24126 


130 ± 12 


13.41 + 0.09 


89.1 + 6.9 


11.5 








1566.428 


0.28853 


111 ± 16 


13.34 + 0.11 


70.7 + 8.5 


11.5 




tentative, S 




1591.389 


0.30906 


78 + 9 


13.20 ± 0.10 


45.3 + 5.0 


11.5 




tentative, S 




1600.816 


0.31682 


62 + 14 


13.05 + 0.17 


43.5 + 6.4 


11.4 




tentative, S 




1604.410 


0.31978 


389 ± 16 


13.99 ± 0.09 


74.4 + 8.7 


12.1 


Ly/?, Ovi 






1610.496 


0.32478 


86 + 12 


13.24 + 0.15 


46.1 ± 10.2 


11.6 




tentative, M 




1676.519 


0.37909 


116 + 19 


13.36 + 0.15 


72.0 ± 8.9 


11.5 




tentative, S 




1680.862 


0.38266 


403 ± 77 


13.92 + 0.42 


200.1 ± 22.8 


11.6 




tentative, M 




1723.654 


0.41786 


88 + 10 


13.25 + 0.08 


50.7 + 3.9 


11.5 




tentative, A 





a la errors are given. 



b Lines listed as tentative are all clearly detected with W r /a > 4.4 except for the z — 0.08587 line toward PG 1116+215 with 
Wr/a — 2.6. However, the tentative lines may not have the required Gaussian line shape to be considered as good candidates 
for single component thermally broadened BLAs. The following abbreviations are used to explain why the detected lines are 
marked tentative: A = line is possibly asymmetric; B = line is blended; M = evidence for multi-component substructure; R = 
possible artifact caused by the MAMA repeller wire; S = S/N in the profile too low to reliably evaluate possible component 
structure. 



P. Richter, B.D. Savage, K.R. Sembach, and T.M. Tripp: Broad Ly a absorbers 



11 



3.5. PG 1116+215 and PG 1259+593 

Detailed analyses of the BLAs towards PG 1116+215 
and PG 1259+593 have been presented by Sembach et 
al. (2004) and Richter et al. (2004), respectively. We here 
briefly summarize and update the findings from these 
studies. Toward PG 1116+215, Sembach et al. (2004) find 
seven BLA candidates along a total, unblocked redshift 
path of Az = 0.133. We regard four of these absorbers as 
tentative detections (see Table 2, last column). With an 
average S/N of ~ 17, the study implies dN^h/dz w 23 
for log (N/b) > 11.2 for the reliably detected systems (~ 
53 including the uncertain cases). Toward PG 1259+593 
there are 29 BLA candidates along relatively long, un- 
blocked redshift path of Az = 0.355 (Richter et al. 2004). 
Following the selection criteria presented above, 22 of 
these candidate BLA systems are tentative detections. 
The average S/N in the the PG 1259+593 spectrum is 12, 
but note that the S/N varies within a factor of ~ 1.5 along 
the spectrum. The number of reliably detected broad 
systems (= 7) corresponds to dN^/dz « 20 for log 
(N/b) > 11.1 — 11.4 (~ 82 including the uncertain cases). 
Column densities, b values and other information about 
these broad Lya candidates towards PG 1116+215 and 
PG 1259+593 are listed in Tables 2 and 3. 



4. Properties of \ow-z broad Ly a absorbers 

4.1. Number distribution 

The measured values for dNBLA/dz along each indi- 
vidual line of sight toward H 1821+643, PG 0953+415, 
PG 1116+215, and PG 1259+593 are summarized in Table 

5. For the definitely detected systems, dN^LA/dz varies 
between 20 (PG 0953+415 and PG 1259+593) and 25 
(H 1821+643). If we include the tentative detections, 
dN B LA/dz ranges between 25 (H 1821+643) and 82 
(PG 1259+593). The discrepancy between the values for 
dNftLA/dz along the individual sight lines may be due 
to intrinsic scatter from the low-number statistics, possi- 
bly reflecting the peculiarities of the WHIM distribution 
along individual sight lines. On the other hand, part of 
the scatter may be caused by the varying sensitivity to 
detect broad Ly a absorption in these spectra due to the 
varying S/N. The PG 1259+593 sight line, for instance, 
provides a large number of tentative BLA detections in 
regions with relatively low S/N. Most likely, a fair frac- 
tion of these tentative BLA candidates are not WHIM 
absorbers but are caused by noise and unresolved compo- 
nent structure. Combining these four lines of sight, we 
have a total redshift path of Az = 0.928 to detected 
BLA candidates. The comoving total path length AX 
sums up to AX = 1.071, as calculated from equation 
(6). Along this redshift path, the STIS data then yield 
a mean value of gL/Vbla / dz = 22 (53) 1 for absorbers with 

1 22 for the secure detections, 53 for entire candidate sample 
including the tentative cases. 



log (N/b) > 11.1 to 11.4, with an absorber weighted aver- 
age of log (N/b) > 11.3. 

4.2. Identification of different absorber types 

We can distinguish between BLAs that are detected solely 
in Lya and systems that show associated absorption in 
Ly/3 and/or Ovi. Associated Ly/3 and/or Ovi absorp- 
tion may be of great importance to gain insight into the 
physical conditions of these systems. In our sample there 
are eight BLA systems that show associated absorption by 
Ly/3 and/or O vi, as indicated in Table 2, seventh column. 
In Table 3 we summarize detailed information about these 
systems, such as Hi and Ovi b values and Ovi velocity 
offsets. 

4.2.1. BLAs with associated Ly/3 absorption 

Five broad Ly a systems in our sample have an associated 
broad Ly (3 line. For the occurrence of associated Ly 3 ab- 
sorption a minimum strength of log (N/b) « 11.8 is re- 
quired. This means that the high-column density systems 
can be detected in Ly/3, if this line falls into a spectral 
region not blended by other ISM/IGM lines. For a ther- 
mally broadened WHIM absorber, the b value measured 
for Ly 3 should be consistent with that of the Ly a absorp- 
tion. If the b values for Ly a and Ly (3 are largely inconsis- 
tent with each other, blending effects or noise features are 
most likely affecting the data and these systems should 
not be considered any further. An example for this is the 
z = 0.30434 system toward PG 1259+593 (see Richter 
et al. 2004, their Fig. 3f), which shows a relatively broad 
Lya line with b w 65 kms" 1 , whereas the Ly/3 absorp- 
tion is narrow with b w 35 kms -1 . For this system we 
can conclude that blending and noise effects but not ther- 
mal line broadening determine the observed large width 
of the Ly a absorption. This system therefore is not in- 
cluded in our sample. A counter-example is the BLA at 
z = 0.31978 along the same line of sight, which has a b 
value of ~ 70 kms -1 consistent with both Lya and Ly/3 
absorption (Richter et al. 2004, their Fig. 3g). Along all 
four sight lines we find nine BLAs with b > 40 kms -1 
that have associated Ly/3 absorption, but only for five of 
these systems is the width of the corresponding Ly (3 line 
consistent with b > 40 kms -1 . Although our statistics are 
poor, this implies that S/N features, blending issues, and 
other processes may contaminate our BLA sample up to 
a level of ~ 50 percent. This value is consistent with our 
estimate from Sect. 2.4, where we analyzed the influence 
of noise in a set of artificial spectra. 

4.2.2. BLAs with associated Ovi absorption 

From the 20 reliably detected BLAs towards H 1821+643, 
PG 0953+415, PG 1116+215, and PG 1259+593, six show 
associated O VI absorption within a velocity range of ±60 
kins" 1 . Only a small fraction of the broad Lya systems is 
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Table 3. BLA systems with associated Ly/3/0vi absorption 



z 


Sight line 


Det. lines 


b T \r n I) T v a 


h vi a 


Avovi 


tj t / u n vt 

^ STL 1 / v_/ V 1 








[kins -1 ] [kms -1 ] 


[kms -1 ] 


[kms -1 ] 




0.22616 


H 1821+643 


HIA1215.7, O vi AA1031.9, 1037.6 


50.8 


15.7 


60 


3.2 




H 1 891 4-fi.dl. 




4^ (S 49 9 


24.1 


q 


1 Q 


0.19360 


PG 0953+415 


HiAA1215.7, 1025.7 


40.3 








0.04125 


PG 1116+215 


HIA1215.7, O vi AA1031.9, 1037.6 


105.0 




35 




0.06244 


PG 1116+215 


HiA1215.7, O vi AA1031.9, 1037.6 


77.3 




20 




0.00229 


PG 1259+593 


HiAA1215.7, 1025.7, OviAA1031.9 


42.1 




48 




0.14852 


PG 1259+593 


HiAA1215.7, 1025.7 


42.1 








0.31978 


PG 1259+593 


HiAA1215.7, 1025.7, Ovi AA1031.9, 1037.6 


74.3 62.7 


19.3 


10 


3.9 



a & values measured only for lines in the STIS wavelength regime. 



expected to have associated O vi absorption, as the occur- 
rence of O vi critically depends on the metallicity of the 
absorbing gas, its temperature, and its total column den- 
sity. Therefore, most BLAs arise in WHIM gas that is too 
hot, to metal-poor, or too diffuse for O vi absorption to oc- 
cur at detectable levels. Some of these "hybrid systems" 
that have good S/N in both Lya and Ovi absorption 
have been analyzed in detail with respect to their physi- 
cal properties, such as temperature, ionization conditions, 
and metallicity (Sembach et al. 2004; Richter et al. 2004; 
Tripp, Savage, & Jenkins 2000). Under the assumption of 
collisional ionization equilibrium, the column densities of 
H i and O vi can be used to constrain the oxygen abun- 
dance in the gas (see, e.g., the z = 0.31978 system towards 
PG 1259+593; Richter et al. 2004). The measured 6 values 
for H i and O Vi also can provide information about the 
influence of non-thermal broadening components in some 
systems. The comparison between broad H I and O vi ab- 
sorption is often difficult, however, as the observed Ovi 
features do not exactly align in velocity space with the 
broad Lya absorption (see also Sect. 2.4). The velocity 
centroids of broad Ly a and O VI in these six systems de- 
viate by 9 — 60 kms -1 (see Table 3). Such velocity offsets 
could be caused by substructures within the filament. It is 
plausible that WHIM filaments are neither perfectly ho- 
mogeneous (in terms of volume densities), nor arc they 
perfectly isothermal (Kawahara et al. 2005). If substruc- 
tures with higher densities exist, the gas may be able to 
cool, creating distinct regions where O vi absorption may 
occur preferentially once the gas reaches the peak temper- 
ature of O vi near 3 x 10 5 K. The radial velocities of these 
denser spots would most likely differ from the center veloc- 
ity of the entire filament, in particular if the gas stays in a 
hydrostatic equilibrium (e.g., Valageas, Schaeffer, & Silk 
2002). Cooler spots then would move towards the mass 
center of the filament, creating a cooling flow. Such effects 
may lead to a non-uniform optical depth distribution of 
H I and O vi along the WHIM filament and could produce 
velocity offsets between the integrated absorption profiles 
of these two species. It therefore remains unclear whether 
the observed broad Hi and the Ovi trace the same gas 
phase. 

Also the comparison between Hi and Ovi 6 values 
provides information about possible non-thermal broad- 



ening effects, as outlined in Sect. 2.3. For three out of the 
six Hi/Ovi absorber pairs accurate Ovi b values can be 
measured and the resulting 6hi/6ovi ratios are listed in 
the last row of Table 3. Only the system at z = 0.31978 
towards PG 1259+593 has 6m ~ 4 6ovi, as expected for 
pure thermal line broadening (equation 2). The systems 
at z = 0.22616 and z = 0.26659 towards H 1821+643 both 
have 6hi < 4 &o vi, thus indicating that 6 non -th dominates 
over 6 t h in these systems. The large widths of the O vi 
lines could be caused by unresolved substructure, while 
the multi-component O vi still could arise in hot gas that 
produces thermally broadened Hi absorption. Since it is 
unclear whether the broad Lya absorption and the Ovi 
arise in the same gas phase, we have to conclude that 
the comparison between BLA and O vi line widths re- 
mains inconclusive. Two of the six broad Lya/ Ovi sys- 
tems occur in the line of sight towards H 1821+643. This 
sight line exhibits a overabundance of intervening O vi 
systems compared to other low-z lines of sight (see Savage 
et al. 2002; Tripp, Savage, & Jenkins 2000). However, we 
here find that the H 1821+643 sight line is not particularly 
rich in broad Lya candidate systems, when compared to 
PG 0953+415, PG 1116+215, and PG 1259+593. 

4.3. CIE properties of the BLAs 

Despite the remaining uncertainty about the interpreta- 
tion of the the observed line widths of BLAs and the un- 
certain contribution of photoionization, we now want to 
estimate the total gas content of these absorbers assum- 
ing the ideal case in which the lines are predominantly 
thermally broadened while the gas is in collisional ion- 
ization equilibrium. We use the measured 6 values and 
H I column densities from our STIS sample together with 
equations (2) and (3) to derive the total hydrogen column 
density, A(H), in these systems. In Table 4 we have listed 
for each of the 49 BLA candidates values for log T, log 
/h, and log iV(H). The BLAs cover a temperature range 
from ~ 10 5 to ~ 2.4 x 10 6 K for pure thermal broadening, 
and the resulting total hydrogen column densities span a 
range from ~ 10 18 to ~ 10 21 cm~ 2 . All of the very broad 
systems with 6 > 100 kms" 1 are considered to be tenta- 
tive detections (see Table 4), so that the validity of the 
highest column density systems is particularly critical. 
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Fig. 6. Number distributions of log 7V(Hl), b, and log 7V(Hn) for our BLA sample are shown. The thick line indicates 
the distribution for the secure cases, the thin line shows the distribution for the entire BLA sample. 



4.4. Distribution of hydrogen column densities and b 
values 

In Fig. 6 we show the distribution of Hi column densities, 
b values, and (estimated) ionized hydrogen column den- 
sities (A(Hii) w Atotai) for our BLA sample. Thin lines 
indicate the distribution for the entire candidate sample 
(49 absorbers), the thick lines indicate the distribution 
for the secure cases (20 absorbers). The Hi column den- 
sity distribution (Fig. 6, left panel) peaks in the range log 
AT(Hi)= 13.2 - 13.4. No system with log iV(Hl)> 14.0 
has been found, demonstrating that the neutral gas con- 
tent of these absorbers is generally low. Note that for log 
N(Hi)< 13.0 our sample is incomplete due to the limited 
sensitivity in the lower S/N regions of the spectra. Only 
a few systems are observed in this range, but it is proba- 
ble that the intrinsic broad Ly a number distribution rises 
further towards the low-column density end. The b value 
distribution of the BLAs (Fig. 6, middle panel) is strongly 
peaked at b = 40 — 60 kms -1 and shows a tail towards 
higher b values. The b value distribution is difficult to in- 
terpret, as there are various selection effects that have to 
be taken into account. First of all, our sample is incom- 
plete for higher b values. Since the sensitivity to detect 
broad Lya lines scales with (N/b) (see equation 4), we 
preferentially detect narrow lines rather than broad lines 
for a given column density level. It is probable that a con- 
siderable fraction of very broad lines lie below our detec- 
tion limit and do not show up in the b value distribution 
shown in Fig. 6. Secondly, if non-thermal line broadening 
mechanisms such as peculiar gas motions are important 
also for the Lya forest lines, we expect to find a number 
of lines from the photoionized Ly a shifted into the b value 
regime between 40 and 100 kms -1 . Such lines may signifi- 
cantly contaminate our broad Ly a sample at low b values, 
whereas the contamination should be less severe at higher 
b values. The number of WHIM absorbers with b values 
b = 40 — 60 kms -1 shown in Fig. 6 thus may be overesti- 
mated. At the sensitivity of the current data, we should be 
able to detect very broad (b > 100 km _1 ) Lya absorption 
toward H 1821+643, and possibly PG 0953+415 down to a 



level of log (N/b) « 11.2 and 11.4, respectively. However, 
all except one of the Ly a lines detected in these two direc- 
tions have b < 100 km -1 . This result stands in stark con- 
trast to the PG 1116+215 and PG 1259+593 sight lines, 
for which several lines with b > 100 km -1 are found. 
Differences in the number of very broad Lya absorbers 
for the four sight lines might be caused by geometrical 
effects resulting from the intersection of the cosmic web 
filaments at different angles. If so, this would imply that 
the line broadening is dominated by flows within the fila- 
ments or by projection effects. Alternatively, the four sight 
lines may intersect regions of the web with very different 
physical conditions, in which case the line- width differ- 
ences may be dominated by the temperature of the gas. 
Distinguishing between these two possibilities will require 
observations of a larger number of sight lines to construct 
a complete line width frequency distribution for absorbers 
in a wide variety of cosmic web environments. 

In Fig. 6, right panel, we show the distribution of the 
ionized hydrogen column density, TV(Hn), as estimated 
for each BLA candidate from equations (2) and (3) (see 
Sect. 4.3 and Tabic 4). If our estimate for A(Hn) is at least 
roughly correct, most of the detected WHIM BLA candi- 
dates have total gas column densities between 10 18 — 10 20 
cm~ 2 . Only three of the securely detected systems have 
column densities larger than 10 19 cm~ 2 , while for the 
more uncertain cases the relative fraction of systems at 
these column densities is significantly higher. This is be- 
cause high A(Hn) systems are expected to have the high- 
est temperatures and thus the lowest neutral gas content. 
As these systems will be particularly broad and shallow, 
they are most difficult to detect. Therefore, our sample is 
incomplete for total gas column densities > 10 19 cm~ 2 . 

4.5. Broad Lya absorbers and structures along the 
line of sight 

Important information about the origin of the broad Ly a 
absorption can be obtained by comparing the redshift dis- 
tribution of BLAs with that of galactic structures along 
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Table 4. CIE properties of broad Lya absorbers 



z 


logT 


log /h 


log JV(H) a 


Status b 


H 1821+643 


0.02642 


5.20 


5.25 


18.49 




0.12221 


5.07 


5.00 


18.23 




0.16352 


5.27 


5.39 


18.58 




0.18047 


5.22 


5.30 


18.51 




0.22616 


5.24 


5.33 


18.80 




0.26659 


5.10 


5.05 


18.70 




PG 0953+415 


0.02108 


5.41 


5.65 


18.77 


tentative 


0.04382 


5.13 


5.11 


18.45 




0.09315 


5.06 


4.98 


18.64 




0.10294 


5.95 


6.55 


19.87 


tentative 


0.17982 


5.10 


5.06 


18.30 




0.19142 


5.06 


4.98 


18.45 


tentative 


0.19360 


4.99 


4.83 


18.74 




PG 1116+215 


0.01635 


5.15 


5.16 


18.55 




0.04125 


5.82 


6.35 


19.59 


tentative 


0.06072 


5.27 


5.38 


18.66 




0.06244 


5.55 


5.91 


19.09 




0.08587 


5.21 


5.28 


18.18 


tentative 


0.09279 


6.03 


6.66 


20.05 


tentative 


0.13370 


5.62 


6.03 


19.30 


tentative 


PG 1259+593 


0.00229 


5.03 


4.91 


18.48 




0.00440 


6.09 


6.74 


20.39 


tentative 


0.00505 


5.29 


5.44 


18.88 


tentative 


0.01347 


5.81 


6.34 


19.58 


tentative 


0.04250 


5.00 


4.84 


17.90 


tentative 


0.05586 


6.10 


6.76 


20.22 


tentative 


0.06810 


5.66 


6.10 


19.24 


tentative 


0.08041 


5.02 


4.90 


17.87 




0.09196 


5.89 


6.45 


19.58 


tentative 


0.10281 


6.37 


7.10 


20.51 


tentative 


0.13351 


5.14 


5.14 


18.15 


tentative 


0.14034 


5.29 


5.42 


18.48 


tentative 


0.14381 


5.91 


6.48 


19.66 


tentative 


0.14852 


5.03 


4.91 


18.82 




0.15136 


5.41 


5.65 


18.97 




0.16647 


5.72 


6.19 


19.53 


tentative 


0.17148 


6.02 


6.64 


19.89 


tentative 


0.17891 


5.77 


6.26 


19.55 


tentative 


0.18550 


5.65 


6.08 


19.26 


tentative 


0.22861 


4.99 


4.83 


18.30 




0.24126 


5.68 


6.12 


19.53 




0.28853 


5.48 


5.78 


19.12 


tentative 


0.30906 


5.09 


5.04 


18.24 


tentative 


0.31682 


5.06 


4.96 


18.01 


tentative 


0.31978 


5.52 


5.85 


19.84 




0.32478 


5.11 


5.07 


18.31 


tentative 


0.37909 


5.49 


5.80 


19.16 


tentative 


0.38266 


6.38 


7.12 


21.04 


tentative 


0.41786 


5.19 


5.23 


18.48 


tentative 



a Total hydrogen column densities are derived under the as- 
sumption that the gas is in collisional ionization equilibrium 
(i.e., additional photoionization is not considered). 
b The tentative status implies that the BLA, although clearly 
detected, may not qualify as a single-component Gaussian 
broadened absorber for the reasons given in Table 2. 



the individual lines of sight. For the lines of sight toward 
H 1821+643, PG 0953+415, and PG 1116+215 redshift in- 
formation about nearby galaxies is available from previous 
studies (Tripp, Lu, and Savage 1998; Savage et al. 2002). 
These galaxy-redshift data have been obtained with the 
3.5m WIYN telescope on Kitt Peak, using the multi-object 
fiber spectrograph (HYDRA). The measurements provide 
redshift information for galaxies in a field of ~ 1 deg 2 size 
around each quasar. Unfortunately, the completeness of 
these galaxy redshift surveys varies substantially among 
the different sight lines, mostly due to the varying stellar 
contamination. However, within ~ 30f0 of each quasar, the 
survey is ~ 40 percent complete for objects with Bj < 20 
(see Tripp, Lu, and Savage 1998). In Fig. 7 we compare 
the redshift distribution of our BLA candidates to the 
distribution of galaxies along the lines of sight towards 
H 1821+643, PG 0953+415, and PG 1116+215 within a 
radius of 30.'0 from the QSO. From a visual inspection 
of Fig. 7 it appears that BLAs are loosely correlated with 
the occurrence of galaxies along the redshift path. One 
or more galaxies are present within a velocity radius of 
3000 kms -1 around each BLA candidate. However, in con- 
trast to the O vi absorbers (e.g., Sembach et al. 2004) the 
BLAs do not appear to be significantly clustered around 
the peaks of the galaxy-redshift distribution. This sug- 
gests that the BLAs also trace only mildly overdense 
WHIM filaments in the vicinity of individual galaxies, 
whereas the O vi systems may be more closely related to 
higher density gas in galaxy groups. However, we may have 
missed broad systems that are part of more complex multi- 
component absorbers, as we have restricted our sample 
to single-component systems. Since multi-component sys- 
tems most likely arise preferentially in the vicinity of in- 
dividual galaxies and galaxy groups due to the more com- 
plex gas distribution in these environments, our BLA sam- 
ple is biased against WHIM absorbers close to galaxies 
and galaxy groups. Therefore, the observed weak correla- 
tion between BLAs galaxies may simply be a result of our 
BLA selection criteria. Additional galaxy-redshift data for 
other lines of sight (e.g., PG 1259+593) will be required 
to improve the statistics and to investigate the relation 
between BLAs and galaxies in more detail. 

4.6. Estimate for f2& 

Following equation (5), the total hydrogen column den- 
sities derived for each of the detected BLA candidates 
toward H 1821+643, PG 0953+415, PG 1116+215, and 
PG 1259+593 (Table 4) allow us to estimate the baryon 
content of the broad Lya absorbers. As discussed ear- 
lier in this paper, there is strong evidence that a signif- 
icant fraction of the detected BLA candidates have ori- 
gins other than the WHIM, and that the b- value distribu- 
tion of the WHIM H I absorbers is affected also by non- 
thermal broadening processes. Moreover, we have shown 
that photoionization from the UV background may sig- 
nificantly influence the ionzation fraction in the BLAs. 
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Fig. 7. Distribution of galaxies within ~ 30.'0 of the di- 
rection of H 1821+643, PG 0953+415, and PG 1116+215. 
The bins have a redshift width of Az = 0.001. The red- 
shifts of the BLAs are indicated by filled circles (secure 
detections) and open circles (tentative detections). The 
smaller circles indicate BLA systems with b < 100 kms" 1 , 
the larger circles indicate very broad systems with b > 100 
kms -1 . 



However, as it is impossible to reliably quantify the con- 
tamination from these alternative processes, we now want 
to provide an estimate for Qb(WHIM) assuming the ideal 
case (i.e., without any contamination and without pho- 
toionization). Summing over the hydrogen column densi- 
ties, log N(R), listed in Table 4, we obtain a total hydro- 



gen column density of log 7V(H)= 20.26 for the reliably 
detected systems, and a value: of log JV(H)= 21.42 includ- 
ing the tentative systems. The total path length to detect 
BLA along the four lines of sight is AX — 1.071 (see 
Sect. 4.1 and Table 5). From equation (6) we then derive 
fi 6 (BLA)= 0.0027 /i7o _1 considering only the 20 reliable 
cases and O b (BLA) = 0.0381 h 70 -1 for all 49 BLA can- 
diatcs. Estimates for Jib (BLA) for each individual sight 
line are listed in the last column of Table 5. Current cos- 
mological models favour a value for of ~ 0.045, and 
the contribution from the WHIM to at z — is ex- 
pected to be ~ 40 percent (e.g., Dave et al. 2001), so 
that fi b (WHIM)« 0.018 is predicted in the local Universe. 
The value for S1;,(BLA) derived from the restricted sample 
therefore is ~ 16 percent of the expected contribution from 
the WHIM to fib- This is slightly higher than the contri- 
bution from intervening O VI absorbers (fif,(0 Vi)> 0.0022 
h 70 - 1 ; Danforth & Shull 2005; Sembach et al. 2004). The 
estimate from the entire BLA candidate sample, however, 
exceeds the expected value for fif,(WHIM) by a factor 
of ~ 2. The reason for this is that the by far dominat- 
ing contribution to fib (BLA) comes from the very broad 
lines, which are exclusively tentative detections. The re- 
liability of these few very broad features is very uncer- 
tain, and the unrealistically high value for fib (BLA) from 
these systems indeed implies that a considerable fraction 
of these features probably are not related to hot gas in the 
WHIM. Note that broad lines with b > 100 kms' 1 make 
up ~ 90 percent of the baryon content in the total sample, 
although they contribute only with ~ 25 percent to the to- 
tal number of absorbers. If we take non-thermal broaden- 
ing processes and additional photoionization into account, 
we believe that our estimate of f2b(BLA)= 0.0027 /i7o -1 
from the limited data sample represents a realistic lower 
limit for the baryon content of broad Ly a systems at low 
z for a detection threshold of log (N/b) « 11.3. 

5. Summary 

In this paper we have studied intervening broad Lya ab- 
sorbers at low redshifts that probably trace large amounts 
of baryonic matter hidden in the warm-hot intergalactic 
medium (WHIM) at temperatures between ~ 10 5 and 
~ 10 6 K. In ionization equilibrium a very small fraction 
(< 10 -5 , typically) of the hydrogen in the WHIM is ex- 
pected to be neutral. Depending on parameters like the 
gas temperature and the total gas column density, a cer- 
tain fraction of this neutral gas should be detectable as 
thermally broadened Lya absorption at column densities 
log A(Hi)< 14 and b > 40 kms -1 . Given the limita- 
tions in sensitivity for detecting very broad and shallow 
absorption features in currently existing STIS data, it be- 
comes clear that only certain regions of the WHIM (those 
with relatively low temperatures and high total gas col- 
umn densities) can be traced with broad Ly a absorption. 
If the gas is in CIE and if thermal broadening dominates, 
one can directly estimate the baryon content of WHIM 
broad Ly a absorbers (BLAs) from the measured H I line 



16 P. Richter, B.D. Savage, K.R. Sembach, and T.M. Tripp: Broad Ly a absorbers 



Table 5. Sight line properties 



Sight line 


4m 


S/N a 


log (iV/6) min b 


AZBLA 


AXbla 


A<BLA 


dN B LA/dz 


n b (BLA) c 


















[xlO- 2 ^,,- 1 ] 


H 1821+643 


0.297 


21 


11.2 


0.238 


0.273 


6(6) d 


25 (25) d 


0.13 (0.13) d 


PG 0953+415 


0.239 


12 


11.4 


0.202 


0.224 


4(7) 


20 (35) 


0.08 (0.55) 


PG 1116+215 


0.177 


17 


11.2 


0.133 


0.144 


3(7) 


23 (53) 


0.22 (2.08) 


PG 1259+593 


0.478 


12 


11.4 


0.355 


0.430 


7(29) 


20 (82) 


0.45 (8.37) 


Total 








0.928 


1.071 


20 (49) 


22 (53) 


0.27(3.81) 



a Typical S/N per resolution clement (~ 7 kms -1 ). 

b Typical detection limit for BLAs. N is in units cm -2 , b in units kms -1 . 

c Estimates for H b (BLA) arc calculated assuming pure thermal broadening for the BLAs and collisional ionization equilibrium. 

d Plain numbers refer to the restricted BLA candidate sample, numbers put in parentheses refer to the total BLA candidate sample (including the 
tentative cases). 



widths and column densities. However, the CIE hypothe- 
sis probably is only a poor approximation for low-density 
absorbers, so that the estimated ionization fractions are 
uncertain by up to 50 percent. In addition, the interpreta- 
tion of BLA line widths is complicated, as non-thermal 
broadening processes like peculiar gas motions, macro- 
scopic turbulence, the Hubble flow, as well as noise fea- 
tures and continuum undulations potentially mimic some 
of the broad absorption features attributed to the WHIM. 
Therefore, every WHIM filament that contains sufficient 
amounts of neutral hydrogen will produce broad Ly a ab- 
sorption, but conversely, not every broad spectral feature 
observed in QSO spectra is related to the WHIM. Due 
to the non-thermal broadening processes, we expect that 
the error rate for detecting the WHIM with broad Lya 
absorption may be as high as 50 percent and that for true 
WHIM BLAs the baryon content may be overestimated 
when relying solely on the observed b value distribution. 
However, it is clear that BLA systems represent an im- 
portant baryon reservoir in the local Universe. Thus, their 
frequency in high-resolution QSO spectra and their phys- 
ical properties should be studied to the greatest extent 
possible. 

Our analysis of the two STIS spectra of H 1821+643 
and PG 0953+415 has resulted in the detection of 14 
broad Lya absorber candidates with 12.88 < log -/V(Hi)< 
13.91 and 40 < b < 122 kms" 1 along a total (un- 
blocked) redshift path of Az — 0.440. These measure- 
ments complement previous analyses of broad Lya ab- 
sorption towards PG 1259+593 (Richter et al. 2004) and 
PG 1116+215 (Sembach et al.2004). Considering all four 
lines of sight, we find 20 reliably detected BLA systems 
and additional 29 tentative cases along a total redshift 
path of Az = 0.928. With 20 (49) detected BLA candi- 
date systems, the number of systems per unit redshift is 
dN^x/dz = 22 (53). Most of these systems have Hi col- 
umn densities between log N = 13.0 — 13.6 and b val- 
ues 40 < b < 100 kms -1 . Assuming CIE and pure ther- 
mal line broadening, the total hydrogen column density in 
these absorbers typically ranges between log N = 18 — 20. 
Unfortunately, most of the systems that have relatively 
high total hydrogen column densities log N > 19 (and 



thus the largest baryon content) belong to the tentative 
detections, and thus are highly uncertain. 

Five of the stronger BLA systems show associated 
broad Ly /3 absorption. The presence of broad Ly f3 is of 
great help to assure that the Hi absorption in these sys- 
tems is indeed intrinsically broad and not caused by noise 
features and continuum undulations. Six BLA candidate 
systems show associated O vi absorption. In all cases, the 
Ovi absorption shows a velocity offset compared to the 
Hi absorption. Moreover, in two out of three cases the 
measured O vi b values are incompatible with the b val- 
ues measured for Hi, possibly implying unresolved compo- 
nent structure and non-thermal line broadening. However, 
as it is not clear whether or not the BLAs and the Ovi 
systems trace the same gas phase, the comparison be- 
tween these two WHIM absorber types remains incon- 
clusive. If we consider only the 20 reliably detected BLA 
candidates, their baryon content sums up to a value of 
n b (BLA)= 0.0027 h^- 1 (assuming CIE), which is - 16 
percent of what is expected to be the contribution of the 
WHIM to Clb at z = and about 30 percent more than 
what is estimated for intervening O Vi systems. If we in- 
clude the tentative detections we derive a high value for 
fib(BLA) of 0.0381 /170 _1 , suggesting that a fair number 
of these tentative detections are not related to hot gas in 
the WHIM but rather are caused by noise features and 
non-thermal line broadening. These findings demonstrate 
that it is quite difficult to exactly pinpoint the baryon 
budget of the BLAs with currently existing STIS data. 
Next to the BLAs that are not related to the WHIM, the 
possible overestimate of fif,(BLA) due to non-thermal line 
broadening and, conversely, the possible underestimate of 
Ob (BLA) due to photoionization effects further compli- 
cate the situation. However, we believe that the value of 
6 (BLA)= 0.0027 /i7o _1 is a reliable lower limit for the 
baryon content of BLAs that lie above a detection limit 
of log (N/b) w 11.3. Our study therefore suggests that 
broad Lya absorbers arising in WHIM filaments contain 
a significant fraction of the baryons at z = 0. 

The BLAs represent an excellent way of studying the 
baryonic content of WHIM gas in the temperature range 
from 10 5 to 10 6 K. Although there are substantial compli- 
cations associated with the separation of the thermal and 
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non-thermal broadening components and the uncertain 
role of photoionization, the fact that the baryon-content 
estimate does not require knowledge of the metallicity of 
the gas is a major advantage for using the BLAs when 
determining the distribution and baryonic content of the 
WHIM. A better understanding of these systems will re- 
quire additional high-resolution ultraviolet spectroscopic 
observations and detailed comparisons with models of the 
hydrodynamical evolution of the gas. Observations of the 
BLAs with the Hubble Space Telescope Cosmic Origins 
Spectrograph (HST/COS) and future ultraviolet spectro- 
scopic facilities may be the only direct way of probing 
low-metallicity regions in the cosmic web. Comparing the 
frequency of observed BLAs with values of dN/dz pre- 
dicted by the simulations would help to constrain the hot 
gas content of the web and test the importance of WHIM 
heat sources other than gravitational collapse of the web 
filaments (i.e., galactic feedback by supernovae). Next to 
the important future X-ray observations of the WHIM 
in emission and absorption, spectral imaging of the Ly a 
emission over large enough fields of view to quantify the 
morphology of the WHIM gas would provide strong tests 
of the formation of large-scale gas structures and their 
relationship to galaxies. The WHIM Lya and metal-line 
emission is expected to be faint but could be observed 
with a moderate cost observatory designed for such an in- 
vestigation (see Sembach et al. 2005b for one such mission 
concept). 
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